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ABSTRACT 



Aims. We consider polarimetric and thermal-emission properties of comet dust and show how and why they can be used 
for classification of comets. 

Methods. We provide a statistical analysis of comet polarimetric, thermal, and orbital characteristics. We perform 
computer simulations of polarization and infrared spectra considering comet particles as ballistic particle-cluster and 
cluster-cluster aggregates (BPCA and BCCA) consisting of submicron spherical grains. 

Results. Comets can be divided into two groups: Type I, characterized by high gas/dust ratio, low polarization, and a 
weak or absent 10 /im silicate feature, and Type II, for which a low gas/dust ratio, high polarization, and strong silicate 
feature are typical. We show that the low polarization is the apparent result of depolarization by gas contamination at 
low dust concentration, which, in turn, results from the dust in Type I comets being concentrated near the nucleus. The 
simulations of thermal emission show that for more porous particles (BCCA), the silicate feature is more pronounced 
than more compact ones (BPCA), for which it even vanishes as the particles become larger. We also show that in both 
types of comets the main contribution to light scattering and emission comes from particles larger than 10 micron. 
Conclusions. The strength of the silicate feature in the cometary infrared spectra suggests that the dust in Type II 
comets consists of high-porosity aggregates, whereas the dust of Type I comets contains low-porosity ones. This is 
consistent with the polarimetric features of these comets, which indicate that the dust in Type I comets tends to 
concentrate near the nucleus. This may result from the predominance of highly processed particles in Type I comets, 
whereas in Type II comets we see pristine or slightly-processed dust. This conclusion is in accordance with the orbital 
characteristics of the comets. We have found that the strength of the silicate feature correlates with the semi-major 
axis of periodic comets and, for short-period comets, with the perihelion distance. Thus, the silicate feature weakens 
due to compaction of aggregate particles if a comet spends more time in the vicinity of the Sun, which allows the comet 
to accumulate a mantle on the surface of its nucleus. 
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1. Introduction W can be estimated directly from photometric or polari- 

metric observations. It characterizes the gas/dust ratio in 
The idea that properties of cometary dust can be a ba- the coma at the moment of the observations and uses the 
sisJot comet classification appeared after |Chernova et al.| same wave i e ngths as those used to measure the polariza- 



Ml) and i Levasseur-Regourd et al.| Um> noticed that ^ whereas rigorous determination of the gas/dust ra- 

comets can be divided into two groups according to the ^ which ig usually defined as the mags ratio of water 

values of the polarization of their dust at phase an- to the dustj requires special infrared observations and has 

gles a « 80 - 100°: the polarization tends to be larger mogt often appeared to be done for different dates from 

than 20% for one group of comets and^maller than 15% when polarization was measured. It turned out that low- 

for the other group. | Chernova et al. | U993D also showed po i arizat ion comets have high values of W, W > 1000 A; 

that high-polarization comets are characterized by low they are namdy charac terized by strong gas emission in 

gas/dust ratios. Quantitatively, they followed the approach comparison with the continuum. This allowed the comets 

of| Krishna Swamy | ^M® and characterized the dust abun- with high w to be called « gas _ rich comets" whereas the 

dance m comets using the ratio of the fluxes in the C 2 other "dust -rich comets" usually have W < 500 A. It was 

band centered at A = 5140 A and the nearby contm- a l so found dChernova et al-lflU ILevasseur-Regourd et al.l 

uum at A = 4845 A. This ratio is denoted as W and |HajmeJl2QQ0) that comet polarization correlates with 



is measured in A, since the flux in the C 2 band is inte- the thei . ma i infrared characteristics: higher polarization is 
grated over the whole band, whereas the flux for the con- usually accompanied by stronger 10 fim silicate feature and 
tinuum is measured for a unit of wavelength. Parameter higher dust temperature 



Send offprint requests to: L. Kolokolova 
email : ludmilla® astro .umd.edu 



In this paper we consider the polarimetric and thermal 
properties of comet dust. We provide their interpretation 
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based on recent observational data and theoretical simula- 
tions of the dust light scattering and emission. We present 
updated data on polarization for dust-rich and gas-rich 
comets and show that they indicate a different distribution 
of the dust within the coma. Our theoretical simulations are 
used to analyze which properties of the dust particles influ- 
ence the polarization and strength of the silicate infrared 
feature. This, together with the orbital characteristics of 
the comets, allows us to reveal the reasons for the existence 
of the two groups (types) of comets. 



2. Polarimetric properties of comet dust: maximum 
polarization 

The dependence of comet polarization on phase angle is 
usually described by (1) position and value of the mini- 
mum (negative) polarization, (2) position of the inversion 
angle (the phase angle where polarization changes from 
negative to positive) and the polarization angular gradi- 
ent, dP/da, at the inversion point, and (3) position and 
value of the maximum (positive) polarization. Amazingly, 
all comets show very similar characteristics of negative po- 
larization (position of the minimum at a = 10° with the 
value P « —1.5%, inversion angle at 20-22°, polarimetric 
slop e dP/da sa 0.3%/deg). Ho wever IChernova et~aTI (jl993l ) 
and lLevasseur-Regourd et all (|l996f ) found that at a > 35°, 
the polarization curve "forks" , forming two groups of curves 
with high and low maximums of polarization. The two 
groups of comets are seen in the second column of Table [T] 
The number s in this column were estimated from the data 
presented in iKiselev et al. I (|2005f ). Very often comets can- 
not be observed at the phase angles of the maximum polar- 
ization, i.e. a « 95°. Therefore in the majority of cases the 
numbers in column 2 represent not the measured but the 
expected values of the maximum polarization, estimated 
from the data obtained at a > 35°. The least reliable data 
are for comets 9P/Tempel 1 and 22P/Kopff, which were ob- 
served at the phase angles whose maximum values reached 
41° and 37° correspondingly. However, at these phase an- 
gles the difference between dust-rich and gas-rich comets 
can already be seen. Both comets showed polarization that 
exceeded the values typical of these gas-rich comets (ap- 
proximately 2-3% higher), and the observed tendency was 
close to the one typical of the dust comets. The table lists 
only those comets whose infrared characteristics are also 
known as we need this for the following analysis (see Sect. 
3). The third column of Table [T] summarizes the data on the 
comet gas/dust ratio as described in the Introduction. The 
fourth column shows the references from which the values of 
parameter W or other characteristics of the gas/dust ratio 
were taken. Table [1] shows that high polarization is typical 
of comets with low values of W, i.e. with a low gas/dust 
ratio. In the following part of this section we explore possi- 
ble reasons for the existence of two groups of comets that 
are different in their polarization and gas/dust ratio. 

We start with analyzing our computer simulations of 
light scat tering by cometary dust. Us i ng the T-matrix tech- 
nique by Mackow ski fe Mishchenkd (|l996f ). we calculated 
light scatte ring by the most real istic type of cometary 
grains (see iKolokolova et al.l 120041 ): aggregates of submi- 
cron particles. We considered random agg regates b u ilt us- 
ing t he ballistic procedure described in iMeakenl (|l983t 
1984). Two types of aggregates, more compact ballistic 



particle-cluster aggregates (BPCA) and more porous bal- 
listic cluster-cluster aggregates (BCCA) were considered. 
For both types of aggregates we managed to obtain the 
correct behavior for the angular and spectral dependences 
of comet brightness and polarization, as well as the cor- 
rect change in spectral char acteristics (color, po larimetric 
color) with the phase angle (jKimura et all 120031 ). We used 
the Halley-type composition of the dust material and the 
constituent particles (monomers) with a radius of 0.1 /jm. 
M ore detailed ca l culati ons and their analysis are presented 
in lKimura et a l. (2006), where we analyzed how the size of 
monomers and their composition influence the angular and 
spectral characteristics of the scattered light. Our calcula- 
tions show (simil ar trends can be also noticed in th e re- 
sults obtained by lLasue fc Levasseur-Regourd I [20061 ) that 
the maximum polarization gets smaller (keeping the correct 
shape of the curve P(a)) under the following conditions: 

— (1) if the real part of the refractive index gets larger; 

— (2) if the imaginary part of the refractive index gets 
smaller in the range within 0.2-0.8, while the real part 
of the refractive index is larger than 1.8; 

— (3) if the monomers in the aggregates get larger. 

Could any of the reasons listed above reasons be responsi- 
ble for the difference in the maximu m polarization in two 
groups of comets? The calculations bv lKimura et al. (2006) 
show that the reasons for changing the maximum polar- 
ization also influence other characteristics of the scattered 
light making them unrealistic. Thus, reasons (1) and (2) 
change the color of the scattered light to blue, whereas the 
color of the comet dust is usually red (blue colors observed 
for some g as-rich comets are a rtifacts from gas contami- 
nation, see lA'Hearn et ailll995f ). Reason (3) leads to some 
incorrect behavior by the negative polarization. We con- 
clude that a difference in comet dust size or composition 
could not explain the difference in comet maximum polar- 
ization without producing a contradiction with the rest of 
the observational facts. Our numerical simulations instead 
suggest that polarization should be intrinsically high re- 
gardless of the properties of the dust particles. 

Recent observational data (IKiselev et al.l 120011 : 
IKiselev et al~l 12004 Uewitt I 12004 Uockers et al.l 120051) 
have shed new light on the two groups of polarimetri- 
cally different comets. They show that low polarization 
manifests itself either when the polarization is averaged 
over a large area in the coma or far from the nucleus. 
Approaching the nucleus, the polarization increases and, 
in the near-nucleus region, reaches the values typical of 
high-polarization comets. The difference in the behavior of 
polarization with the distance from the nucleus for comets 
with different polarimetric properties is seen in Fig. [1] As 
Table [1] demonstrates, comets with low polarization are 
also characterized by the p redominance o f gas over dust 
in the coma. Based on this. IKiselev et ah! (|200ll ) and then 
iJockers et al.l (|2005| ) explain the low average polarization 
by the contribution of gas emission in the continuum, 
which causes the depolarization of the continuum light. 
Their analysis was based on a detailed comparison of 
the wavebands of narrow-band filters and corresponding 
parts of almost-simultaneously-obtained cometary spectra, 
and it allowed direct determination of gas contribution 
into continuum. They find that as we get closer to the 
nucleus, concentration of the dust increases. This decreases 
the relative contamination by gas emission and makes 
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Table 1. Comet gas/dust characteristics and polarization. 



Comet name 


Polarization^ 


W 


Reference on 




at a « 95°~ 


A 


gas/dust ratio 


2P/Encke 


~ 8 - 10%" 


weak continuum 


A'Hearn et al. (1995) 


C/1975 Nl (Kobayashi-Berger-Milon) 


ra 10% " 


weak continuum 


Winiarski et al. (19921 


23P/Brorsen-Metcalf 


10-15% 


1840-2010 


Chernova et al. (19931 


27P/Crommelin 


10-15% 


860-1470 


Chernova et al. ( 19931 


C/1989 Ql (Okazaki-Levy-Rudenko) 


« 16% 


gas- rich comet 


Winiarski et al. (19921 


D/1996 Ql (Tabur) 


12-17%^ 


gas-rich comet 


Kawakita et al. ( 19971 


C/1989 XI (Austin) 


< 20% 


180-1300 


Joshi et al. (1992) 


67P / Churyumov-Gerasimenko 


< 20% 


relatively dust-rich comet 


Hanner et al. (1985b) 


9P/Tempel 1 


«24% 


dust-rich comet 


KiiDDers et al. (20051 


2 IP/ Giacobini-Zinner 


« 24% 


50 


Chernova et al. (19931 


C/1983 HI (IRAS-Araki-Alcock) 


>25% d 


dust-poor comet 


Hanner et al. (1985a1 


22P/KopfI 


» 26% 


30 


Chernova et al. (19931 


lP/Halley 


» 26% 


100 


Chernova et al. (19931 


C/1987 PI (Bradfield) 


« 26% 


160 


Chernova et al. (19931 


C/1990 Kl (Levy) 


« 26% 


220 


Chernova et al. (19931 


C/1996 B2 (Hyakutake) 


« 26% 


200 


Kiselev & Velichko (19981 


C/1975 VI (West) 


» 26% 


strong continuum 


Rosenbush (19851 


C/1995 01 (Hale-Bopp) 


> 30% 


120-160 


Kiselev fc Velichko (19971 



a the data for the red spectral range 

b as observed through wide-band filters and large diaphragm ([Jockers et al.l l2005l 

c at large diaphragms ( Jock ers~lll997l 1 

d at small diaphragm 



this difference is a property of the comet itself (and may 
indicate some sub-classes of comets), independent of the 
aperture or filter passband. In the case of gas-rich comets, 
the polarization depends significantly on the aperture 
of the measurements and filter passband (Fig. [21 right 
panel), thus these comets cannot be characterized by a 
single "typical" value of the maximum polarization. 

Morphological features such as jets and shells can show 
somewhat different polarization from the polarization of the 
rest of the coma. However, their polarization does not con- 
tribute significantly to the integral polarization of comets, 
especially at large phase angles. This can be seen from the 
fact that the dependences of polarization vs. phase angle 
for dust-rich comets are stable and do not demonstrate any 
noticeable scattering of the data. The more noticeable scat- 
tering of the data for gas-rich comets can be explained by 
the different contribution of gas for dif ferent comets (see, 
iKiselev et all 12004 Uockers et al.ll2005l ). In summary, the 
difference in the comet polarization results not from the 
different size or composition of the dust grains, but from 
the fact that in low-polarization comets the dust is concen- 
trated near the nucleus, leaving distant parts of the coma 
dominated by gas, which contaminates the continuum and, 
thus, depolarizes it. However, this explanation does not an- 
swer why comets form two groups. It just shifts the ques- 
tion about two polarimetric groups of comets to another 
one: why are there two groups of comets that differ in the 
distribution of the dust in the coma? We try to find the 
answer in the following section using computer simulations 
of the thermal properties of cometary dust. 

3. Infrared properties of comet dust: strength of 
the silicate feature 

As mentioned in the Introduction, the values of maximum 
polarization correlate with some thermal characteristics 
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Fig. 1. Comet polarization in the red spectral range versus 
the distance from the nucleus for dust-rich (filled symbols) 
and gas-rich (open symbols) comets. The polarization is 
averaged over the aperture of radius p and is normalized to 
the value obtained at the smallest aperture. 



the polarization values similar to those observed for the 
comets with dust-dominated comae. The coma of gas-rich 
comets gets dominated by gas at nucleocentric distances 
> 1000 km; this is why at the observat ions with low spatial 
resolution (see, e.g. lKellev et al. II2004D the sp at ial gr adient 
in polarizatio n reported by I Jockers et alJ ([2005) and 
iJewitt I (|2004|) may not be noticed. If we make a plot that 
includes the data for different distances from the nucleus, 
the comets still can be divided in two groups. However, 
the difference between the groups will be in the different 
distributions of the polarization over the coma. In dust-rich 
comets, polarization does not show significant changes 
with the distance from the nucleus. Thus, dust-rich comets 
can be characterized by some typical value of maximum 
polarization (Fig. [21 left panel) . Although some difference 
between comets can be noticed, e.g. higher polarization 
for Hale-Bopp or lower polarization for Giacobini-Zinner, 
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Table 3. Complex refractive indices (m 
wavelength A for the synthetic materials. 



i k) at a 



20 40 60 80 100 120 
Phase angle, deg 



40 60 80 100 120 
Phase angle, deg 



Fig. 2. Phase-angle dependence of polarization for dust- 
rich (left panel) and gas -rich(right panel) co mets in the red 
spectral range (based on lKiselev et al. II2005D . The solid line 
shows the best fit curve for dust-rich comets. At large phase 
angles, polarization averaged over a large aperture for gas- 
rich comets (large symbols) is substantially smaller than 
for dust-rich comets illustrating the original idea about two 
groups of comets. Polarization of gas-rich comets measured 
in the near-nucleus area of the coma and properly reduced 
for gas contamination is significantly larger (small sym- 
bols). Vertical lines show the range of polarization change 
as the distance from the nucleus (aperture) changes. 



of comets, e.g. dust temperature and the strength of the 
10 /xm silicate feature. Since there is a strong correlation be- 
tween the temperature and strength of the silicat e feature 
(see, e.g.. iGenrz fc Nev I Il992l; iLisse et al~ll2002ft . we will 
characterize comet thermal emission only by the strength 
of the silicate feature, considering the dust temper ature 
as a redund a nt ch aracteristic. Following iLissd (|2002f) and 
ISitko et al. I {2004), we define the strength of the silicate 
feature as the ratio of the flux between 10 and 11 micron 
to that of the underlying continuum. Although the infrared 
feature of crystalline silicate differs from that of amorphous 
silicate, this does not influence the strength of the silicate 
feature; a detailed analysis of the silicate mineralogy is the 
subject of our separate study. 

Table [5] summarizes the data on the strength of the sil- 
icate feature and demonstrates that comets can be divided 
into two groups, which we call Type I (the comets with a 
weak silicate feature) and Type I I (the comets with a s trong 
silicate feature) comets following [Genrz fc Nev I (119921) . The 
first idea that comes to the mind is that the strong sil- 
icate feature indicates high abundance of silicates in the 
dust, whereas the low silicate feature shows their short- 
age. However, comparing the second columns in Table [1] 
and Table one can see a correlation between the values 
of maximum polarization and the strength of the silicate 
feature: Type I comets are also characterized by low polar- 
ization whereas Type II comets have high polarization. In 
Sect. 2 we showed that polarimetric properties do not in- 
dicate any significant compositional difference between the 
dust in different comets. 

To find the reason for the difference in the infrared 
properties of the comet dust, we calculated infrared spec- 
tra of aggregated particles choosing their size, composi- 
tion, and structure the same as provided the best fit to 
the observed brig htness and pol a rizati on in the calcula- 
tions presented in iKimura et alj |2006), which were dis- 
cussed in Sect. 2. We built our particles as ballistic ag- 
gregates and, to account for different porosity, considered 



Wavelength 

M 



Complex refractive index 
m = n + i k 



8.0 


2 67 4- i 6 1 1 x 10" 


i 


8.5 


2.70 + i 6.23 x 10" 


l 


9.0 


2.71 +i 6.33 x 10" 


i 


9.5 


2.70 + i 6.58 x 10" 


i 


10.0 


2.71 +i 7.82 x 10" 


i 


10.5 


2.84 + i 7.12 x 10" 


i 


11.0 


2.72 + i 7.61 x 10" 


i 


11.5 


2.94 + i 1.16 x 10" 


l 


12.0 


3.16 + i7.71 x 10" 


i 


12.5 


3.14 + i6.07 x 10" 


i 


13.0 


3.03 + i 6.08 x 10" 


i 



BCCA and BPCA particles. Calculations were performed 
for an aggregate made of spherical monomers of radius 0.1 
fim whose composition corresponds to the Halley-comet 
material, i.e. we assume that cometary dust is composed 
of magnesium-rich olivine, organic refractory, amorphous 
carbon, an d pyrrhotite with th eir volume fractions de- 
rived from IKimura et alj (|2006f) . We too k the refractive 
indice s for magnesium-rich olivi ne from | Mukai fc Koikd 
(fl990h - organic refractor y from iLi fe Greenberel (Il997l ). 



1991D, and 



amorphous car bon from iRouleau fc Martini 
pyrrhotite from lBegemann et alJ (1993), with the refractive 
indices of pyrrhotite bel ow A = 10 /im extrapola ted from 
the near- infrared data of lEgan fc Hilgemanl (|l977f ). Tabled 
presents the effective refractive indices of this synthetic 
mixture as a function of wavel ength calculated using th e 
Maxwell Garnett mixing rule (iBohren fc Huffman! Il983h . 
We estimate the absorption cross sections for model aggre- 
gates in the infrared wavelengths using Mie theory along 
with the Maxwell Garnett mixing rule (see iMukai et al.l 
fl99l) . Thus, we calculate thermal emission for a spherical 
particle with the radius chosen to provide a characteristic 
size of the corresponding aggregate. The refractive index of 
the sphere m aterial is calculated using t he Maxwell Garnett 
mixing rule (jBohren fc Huffmanlll983l ) when the material 
was presented as a mixture of the aggregate material from 
Tableland the voids, so that the ratio of the volume occu- 
pied by the material to the volume of the voids corresponds 
to those in the aggregate. 

Our numerous tests showed that the spectral depen- 
dence of the thermal emission calculated with the Maxwell 
Garnett approximation keeps all the features of the spec- 
tral dependence obtained with more rigorous methods. At 
the number of monomers N = 2 10 , an error in the infrared 
absorption cross section arising from this approximation is 
evaluated using the super position T-matrix method, which 
gives a rigorous solution (jMackowski fc Mishchenkol ll996). 
For an alternative error estimate, we also calculate the in- 
frared absorption cross sections for aggregates consisting of 
N = 2 15 monomers us ing the discrete dipole approximation 
(|Okamoto et al.|[l994[) . We find that errors are nearly inde- 
pendent of wavelength, and thus our qualitative discussion 
on the infrared spectra can be validated. 

Figure [3] presents the absorption cross sections divided 
by the volume of aggregate particles composed of the syn- 
thetic mixture of magnesium-rich olivine, organic refrac- 
tory, amorphous carbon, and pyrrhotite. Each panel shows 
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Table 2. Infrared and orbital characteristics of comets. 
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Comet name 


Silicate 


Perihelion, 


Eccentricity 


Semi-major 


Reference 




feature strength^ 


AU 




axis, AU 


for infrared data 


2P/Encke 


<1.1 


0.33 


0.85 


2.2 


Lissc (2002) 












Gehrz et al. (1989) 


/"I /-I r\i—rf AT ~t /TZ 1 1 " T"1 TV X'l 

C/1975 Nl (Kobayashi-Berger-Milon) 


1.0 


0.425 


1.00097 


CO 


AT / -t A O A \ 

Nev (1982) 


23P/Brorsen-Metcalf 


1.0 


0.478 


0.97 


15.9 


Lvnch et al. ( 1992) 


27P/Crommelm 


1.1 


0.743 


0.92 


9.29 


Eaton & Zarnccki (1985) 


C/1989 Ql (Okazaki-Levy-Rudenko) 


1.2 


0.642 


1.0000197 


CO 


Sitko et al. (2004) 


D/1996 Ql (Tabur) 


1.3 


0.842 


0.999474 


1600 


Harker et al. (1999) 


C/1989 XI (Austin) 


1.13 


0.349 


1.000225 


oo 


Sitko et al. (2004) 


67P/Churyumov-Gerasimenko 


1.0 


1.292 


0.63 


3.49 


Hanner et al. (1985b) 


9P/Tempel 1 


<1.1 


1.5 


0.52 


3.13 


Lissc (2002) 


C21P/Giacobmi-Zinner 


1.1 


1.034 


0.71 


3.56 


Hanner et al. (1992) 


/1 AOQ TJ1 /T"D AG A „„1 ' A 

U/1983 xlt (IKAb-Araki-Alcock) 


i n a 
1.1)4 


a aai 
0.991 


u.yyoiio 


1 AA O 

1UU.2 


bitko et al. ( 21)04 ) 


22P/Kopff 


1.2 


1.58 


0.54 


3.43 


Lisse (2002) 


lP/Halley 


1.6 


0.587 


0.97 


19.6 


Harker et al. ( 1999) 


C/1987 PI (Bradfield) 


1.9 


0.869 


0.999697 


2868 


Harker et al. (1999) 


C/1990 Kl (Levy) 


1.8 


0.939 


1.000417 


oo 


Harker et al. (1999) 


C/1996 B2 (Hyakutake) 


>1.5 


0.230 


0.999758 


950.4 


Lisse (2002) 


C/1975 VI (West) 


2.0 


0.196 


0.999971 


6752 


Nev (1982) 


C/1995 Ol (Hale-Bopp) 


2.16 


0.914 


0.995069 


185.4 


Sitko et al. (2004) 


Coefficient of correlation between strength of the silicate feature and corresponding orbital characteristic 


All comets 




-0.193 








Periodic comets 




-0.352 




0.662 




Short-period comets 




0.718 




0.751 





a Although no distinct dependence of the strength of the silicate feature on the heliocentric distance (Harker et al. |[l999l ; 
ISitko et al. 1120041 ) was o bserved, in the case where several values were available we selected the value at the heliocentric distance 
~ 1 AU. Note also that lGehrz et al~l (2005) find no difference between pre- and post-perihelion silicate features. 
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Sphere 
BPCA 
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N=2" 




N=2' 




N=2' 




N=2~ 




N=2^ 




N=2 J 




8 10 12 8 10 12 8 10 12 8 10 12 8 10 12 8 10 12 8 10 12 

Wavelength (jim) 

Fig. 3. Absorption cross section per volume as a function of wavelength for aggregate particles composed of a synthetic 
mixture of magnesium-rich olivine, amorphous carbon, and pyrrhotite. The number of monomers increases from N = 2° 
(left) to 2 30 (right), which correspond to the radii of aggregates ay — 0.1-100 /im. Solid lines: results for particles grown 
under the ballistic cluster-cluster aggregation (BCCA) process; dotted lines: results for particles grown under the ballistic 
particle-cluster aggregation (BPCA) process; dash-dotted lines: results for compact homogeneous spherical particles. 
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Table 4. The porosity P of BPCA and BCCA at various 
sizes. 





H 


act 


[//ml 


1 


- P c 


N 


BPOA 


BOCA 


BPOA 


B(J(JA 


2 U 


0.10 


0.10 


0.10 


1.0 


1.0 


2 5 


0.32 


0.59 


0.76 


0.15 


0.072 


2 io 


1.0 


1.9 


4.6 


0.15 


0.011 


2 15 


3.2 


6.1 


28.0 


0.15 


1.6 x 10 -3 


2 20 


10.0 


19.0 


170 


0.14 


2.3 x 10" 4 




32.0 


62.0 


1000 


0.14 


3.3 x 10~ 5 


2 30 


100 


200 


6000 


0.14 


4.9 x 10~ 6 



a The radius of volume-equivalent spheres. 
b The characteristic radius of aggregate particles: ac = 
/5/3a g where a g denotes the gyration radius of the aggregates. 
c 1 — P = N (a/ac) 3 where a is the radius of monomers. 



the results for aggregates of different sizes, defined by the 
number of monomers N. One can see that the infrared spec- 
tra for those aggregates smaller than N < 2 20 (ay < 10 /im 
where ov = aiV 1 '' 3 ) always show the silicate features irre- 
spective of their sizes and porosities. When the aggregates 
reach radii ay > 10 ^m, low-porosity (Table S]) BPCA par- 
ticles do not show the silicate features, but high-porosity 

(Tabic 2]) BCCA particles retain the f eatures. 

In. the previous studies (see, e.; a.. lHannerl [2fM iLissd 

I2002t iLisse et al. Il2002t iLisse etal.1 12004ft . weak or absent 
silicate features and corresponding low polarization were 
explained by the predominance of large particles in the 
coma, whereas strong silicate features and high polariza- 
tion were presumed, indicating small particles. This con- 
clusion might be correct if the low polarization in comets 
with a weak silicate feature were an intrinsic property of the 
dust. However, as we showed in Sect. 2, all comets in the 
near-nucleus area have similar (high) values of polarization. 
This could be im possible if the particles h ave different sizes 
(see Sect. 2 and lKimura et al.ll2003l.l2006ft . As we concluded 
earlier, two types of comets are characterized by different 
distributions of the dust in the coma: the dust tends to 
concentrate near the nucleus for low-polarization comets. 

Thus, we have to look for another characteristic than 
size of the dust particles that can be responsible for the dif- 
ferent strength of the silicate feature. The analysis of Fig. [3] 
allows us to suggest that the porosity of the dust particles 
may be the very crucial parameter that divides comets into 
two groups. Dust particles are porous in the comets that 
demonstrate a strong silicate feature and more compact in 
the comets with a weak or absent silicate feature. This con- 
clusion is consistent with the different spatial distribution 
of the dust in the comets of different types. Indeed, com- 
pact large particles in the "low-polarization" comets are 
heavy and cannot be accelerated by the gas drag, whereas 
porous, albeit large, particles in "high-polarization" comets 
can easily be moved out of the nucleus. 



4. Physical properties of comet dust: typical sizes 

Before reasoning out the difference in dust porosities, we 
look into a typical size of aggregates. We have noticed from 
our numerical calculations of linear polarization that large 
aggregates are required to produc e the observed negative 
polarization at small phase angles (|Kimura fc Mann Il200l 
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Fig. 4. Computed minimum polarization as a function of 
the size of aggregate. 



iKimura et aLll2003l 12006ft . In Fig.H we summarize our pre- 
vious results for the dependence of the minimum polariza- 
tion on the size of aggregates at a wavelength of 0.6 /zm, 
including our unpublished results. This clearly shows that 
the observed negative polarization of approximately —1.5% 
cannot be attained by small aggregates of a few microns or 
less. To be consistent with both polarimetric and infrared 
observations of comets, we have to admit that the dust 
particles, which provide the major contribution in the light- 
scattering and thermal-emission characteristics of comets, 
are larger than ay — 10 /zm for any comet. This conclu- 
sion is also consistent with the in situ data for the mass 
distri bution of the dust in comet Halley ( McDon nell et al71 
Il987l ). Comet Halley is a typical high-pol arizat i on, str ong- 
silicate-feature comet, i.e., according to iLissd (|2002ft . its 
dust should be characterized by predominantly small parti- 
cles. However, the size distribution of cometary dust derived 
from in situ data on comet lP/Halley shows that large par- 
ticles contribute the maj o rity of the total cro ss-section area 
(|McDonnell et al. I Il987t iFulle et al. I l2000h . We updated 
the distribution of the cross-section area of the particles 
(the parameter, which defines the effectiveness of the light 
scattering and emission) shown in lMcDonnell et al. I (|1987j) 
considering comet particles as BCCA and BPCA with the 
density of the monomer material equal to 2.4 x 10 3 kg m~ 3 
(based on the composition described in the beginning of 
Sect. 3). The results are shown in Fig. [5] One can see 
that the cross-section distribution for BPCA has a local 
maximum around ps 10~ 13 kg (ay ps 3 /im), but reveals 
the importance of much larger particles having radii of 
ay > 10 3 /im. For BCCA particles, the cross section is 
smoothly increasing with the size of particles, indicating 
that large particles make a much more significant contribu- 
tion into light scattering and emission. 
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Fig. 5. Size distr ibution from the in si tu measurements 
for comet Halley (McDonn ell et al71l!987f ) and correspond- 
ing cross-section area of the dust particles for BPCA and 
BCCA. 



5. Synthesis of the observational data and orbital 
characteristics of comets 



In this section we try to find out why the dust in some 
comets is dominated by compact particles, while the par- 
ticles are porous in other comets. Following the way sug- 
gested by numerous previous studies, we check if this comes 
from the in different evolution of the comets, i.e. if there 
are any dynamical characteristics of the comets that corre- 
late with the properties of their dust. lA'Hearn et a l. (1995) 
noticed that the perihelion distance co rrelat es with the 
;as/d ust ratio. According to iLissd (|2002f ) and iLisse et al. I 
2002), there is a correlation between the value of perihe- 
lion distance and strength of the silicate feature. This is 
why perihelion distance (third column in Table [5]) is the 
orbital characteristic whose correlation with the strength 
of the silicate feature we are checking first. Table [3] also 
shows eccentricity, e, of the considered comets. This al- 
lows us to distinguish between new (e > 1), long-period 
(0.99 < e < 1), Halley-type (0.99 < e < 0.9), and short- 
period (e < 0.9) comets. 

The bottom of Table [5] shows the coefficient of correla- 
tion between infrared and orbital parameters. One can see 
that the correlation between perihelion and strength of the 
silicate feature is rather significant for short-period comets 
but loses its significance if long-period comets are included 
in the analysis. Based on this we suppose that the closeness 
of a comet to the Sun is a key issue; thus, not the perihe- 
lion, but the comet average distance from the Sun (or the 
semi-major axis of its orbit) may be the characteristic to 
consider. The fifth column of Table O shows the values of 
this characteristic and its correlation with the strength of 
the silicate feature. This is also illustrated by Fig. [5] We can 
see that all periodic comets show a significant correlation 
with this characteristic. For short-period comets the cor- 
relation with the semi-major axis is even more significant 
than the correlation with the perihelion distance. 

The results of the performed statistical analysis, to- 
gether with the conclusions of Sect. 3, are summarized in 
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Fig. 6. Dependence of the strength of the silicate feature 
on the semi-major axis of the comet orbit for Jupiter fam- 
ily (o), Halley-type (+), and long-period (x) comets. The 
comets from Table 2 are shown by larger symbols. 

Table 5. Classification of comets based on their polarimet- 
ric, infrared, and orbital characteristics 



Characteristic 


Type I 


Type II 


Gas/dust ratio 


high 


low 


Polarization 






- averaged over 


low 


high 


large aperture 






- spatial distribution 


increases as 


rather 




approaching 


homogeneous 




the nucleus 




Silicate feature strength 


low 


high 


Average distance 


small 


large 


from the SunJ^ 






Specifics of dust particles 


compact 


porous 



based on comet dynamical type and value if its semi-major 



Table They allow us to suggest that the more time a 
comet spends in the vicinity of the Sun and the closer it 
approaches the Sun, the more compact are its dust parti- 
cles. 

It is well known that a dust mantle is formed on the sur- 
faces of cometary n uclei along with rele ase of gas and dust 
near the Sun (e.g., iHoupis et al.lH985l ). At the formation 
of the dust mantle, sublimation of volatiles in cometary 
nuclei may lower the porosity of dust aggregates by the 
packing effect, whic h results from anisotro pic sublimation 
of the volatiles (cf. iMukai k Fechtieiri983D . Alternatively, 
high-porosity aggregates, such as BCCA particles, may be 
selectively removed from a cometary nucleus owing to their 
higher mobility in a gas flow compared to low-porosity ag- 
gregates, such as BPCA particles, which tend to stay in the 
nucleus. This produces a deficit of high-porosity aggregates 
in cometary nuclei, which increases as the influence of the 
solar radiation gets stronger and/or longer. 

Because the growth of dust mantles on the surfaces of 
comets depends on their orbital parameters, short-period 
comets are likely to have low-porosity aggregates on av- 
erage, while long-period comets retain high-porosity ag- 
gr egates. This not on ly explains the fact, noticed earlier 
by ISitko et al. I (|2004l ). that most long-period comets show 
silicate features in the mid-infrared spectra, while the in- 
frared spectra of short-period comets are almost feature- 
less; but it is also consistent with correlation between the 
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strength of the infrared feature (as well as polarization 
and gas/dust ratio) and the semi-major axis for all peri- 
odic comets and similar correlation with the perihelion dis- 
tance for short-period comets. This is also consistent with 
the strong correlation be tween the strength of the silicate 
feature and jet activity (jHadamcik fc Levasseur-Regourd I 
2003: lGehrz et al71l2005T ) . since jets from active areas indi- 
cate ejection of the material from the inside, beneath the 
surface mantle. 

Note that new comets, which we could not include in 
our analysis because their semi-major axes cannot be de- 
fined, show a variety of values for the strength of the silicate 
feature (TableEJ. This may indicate some initial features of 
their parent body, which, probably, become less noticeable 
after several close approaches to the Sun. 
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